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Abstract 

We study approximation error bounds of isogeometric function spaces 
on a specific type of singularly parameterized domains. In this context 
an isogeometric function is the composition of a piecewise rational func¬ 
tion with the inverse of a piecewise rational geometry parameterization. 
We consider domains where one edge of the parameter domain is mapped 
onto one point in physical space. To be more precise, in our configuration 
the singular patch is derived from a reparameterization of a regular trian¬ 
gular patch. On such a domain one can define an isogeometric function 
space fulHlling certain regularity criteria that guarantee optimal conver¬ 
gence. The main contribution of this paper is to prove approximation error 
bounds for the previously defined class of isogeometric discretizations. 


1 Introduction 


Isogeometric analysis as presented by Hughes et al. (2005) starts from a dis¬ 
cretization based on a B-spline or NURBS geometry parameterization. The 
standard theory, as developed in Bazilevs et al. (2006) relies on the regularity of 
the geometry mapping. If the geometry mapping is singular, the standard the¬ 
ory does not apply. Therefore, in the standard setting, the underlying geometry 
has to be diffeomorphic to a rectangle. Singular parameterizations have been 
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studied in the context of isogeometric analysis in e.g. Lu (2009); Cohen et al. 


(2010); Lipton et al. (2010) The regularity theory has already been partially ex¬ 
tended to singularly parameterized domains in Takacs and Jiittler (2011, 2012). 
However, the study of approximation properties of isogeometric spaces over sin¬ 
gular patches still remains an open problem, which we study in more detail in 
this paper. In the context of finite elements there exist several related studies 


concerning degenerate isoparametric elements, such as Jamet (1977); Acosta 
and Monz6n| (|2006[). 


In Section we present the notation and the underlying setting that we 
consider throughout this paper. In Section we introduce the hierarchical 
mesh as well as the corresponding function space on it. In Section 3.3 we 


define an L^-stable projection operator onto the hierarchical function space. 
We briefly mention how the construction can be generalized to locally quasi¬ 
uniform refinements in Section |3.4| and present approximation error bounds in 
Section]^ Finally, we conclude the paper in Section 


2 Preliminaries 

Isogeometric function spaces V, as they are present in isogeometric analysis 


introduced in Hughes et al. (2005), are built from an underlying B-spline or 


NURBS space. Hence, to introduce the notation needed, we start this prelimi¬ 
nary section with recalling the notion of B-splines and NURBS. We do not give 
detailed definitions here but refer to standard literature for further reading, see 


Piegl and Tiller 

(1995) 

Prautzsch et al. 

(2002 


2.1 B-splines, NURBS and isogeometric functions 

Throughout this paper we consider the simple configuration of an isogeometric 
function space based on a single NURBS patch with uniform knots. Note that 


this simplification is not necessary, as we point out in Section 3.4 


For simplicity we consider uniform B-splines over the parameter domain 
]0,1[. Given a degree p € and a mesh size h = 1/2", with n S Zj, the Uth 
B-spline, for z = I,..., 2" -|- p, is denoted by &/(s). Here n refers to the level of 
(dyadic) refinement. We assume that the knot vector is open and has uniformly 
distributed interior knots. We moreover denote by 5"^ (s,t) the product of 6/(s) 
and bj(t). Let be the tensor-product B-spline space spanned by b'^j{s,t). 
We have given a weight function w G Sn* on some coarse level n*, with w > 0 
uniformly. For all n> n* we define the NURBS space as A//, = {— : Vn G Sn}■ 
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Moreover, let G S be a geometry parameterization, with 

G : B = ]0,1[^ ^ C 

For n>n*, the h-dependent spaces V„ of isogeometric functions over the open 
domain 17 = G(B) are defined by 

V„ = {ipn : 17 K I </?„ = /„ o G”\ with /„ e A/"™} . 

The goal of this paper is to prove approximation error bounds for h-refined 
isogeometric functions over singularly parameterized domains. In the following 
section we introduce a simple configuration of singularly parameterized domains. 


2.2 Singular tensor-product patches derived from trian¬ 
gular patches 


In this section we construct smooth isogeometric function spaces on singular 
patches, where the singular patches are derived from triangular Bezier patches. 
The configuration presented here is developed in more detail in Takacs (2015), 
which is based on Hu (20011. 

Let u be the mapping 


u: ]0,lp 
{s,tf 


A = {{u,v) : 0 < u < 1, 0 < u < u} 
(s,s7)^ , 


and let F : A —)■ be a regular mapping. We assume that the parameterization 

G is given as 

G = F o u. 

All mappings u, F and G are defined on an open parameter domain and can 
be extended continuously to the boundary. Here, G is singular at a part of 
the boundary, i.e. det VG(s,7) = 0 for all (s,i) S {0} x [0,1]. Let be the 
isogeometric B-spline space on A, i.e. 

= {g}n : A K I (/?„ = /„ o with /„ e 5„} , 


where the inverse of the singular mapping u is equal to 

u-^{u,v) = . 

The various introduced mappings and domains are depicted in Figure 
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Figure 1: Mappings F, u and G for an example domain 


2.3 Regularity conditions 

To prove approximation properties on il, we need the following. 

Assumption 1. There exists a constant Cp, depending only on F and on the 
degree p, such that 

^ ll‘^ o 

for all 0 < k < p + 1 and for all ip G 

Considering a mapping that is continuous of order p on the closure of the 
open domain, we have the following. 

Proposition 1. If F G and detVF > c > 0, then Assumption^ is 

fulfilled. 


The proof of this proposition can be carried out similarly to the proof of 
Lemma 3.5 in Bazilevs et al. (20061. Here, the space is defined as follows. 


Definition 1. Given an open domain D. The space of ^^^-continuous 

functions on the closure of D is defined as the space of functions (^ : I? —>■ K 
with ip G C^{D) such that there exists a unique limit 


y™ dxf^dxf-^ 

yeD 1 2 


dxf^dxf^ 


for all X G dD — D\D and for all |a| = 01+02 < k. Here C^{D) is the standard 
space of fc-times continuously differentiable functions on the open domain D. 
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3 Spaces and operators on the triangle 


In this section we introduce a hierarchical mesh refinement and corresponding 
spline spaces >V„ on the triangle A. In this configuration >V„ is a subspace of 
yV„. Then we introduce a projection operator 11^ : T^(A) —>• >V„, that is 

locally bounded in 


3.1 Hierarchical mesh and refinement 


Let Tn be the Bezier mesh corresponding to the function space >V„ on A, i.e. 
it is a regular mesh mapped by u. We consider a sequence of meshes T„, such 
that To = {A} and 


where 


— {(5 : d S Tn n (A\Ai/2 )} U {d/2 : d S Tn_i}, 


Aj = {(u, v) : 0 < u < j, 0 < V < u} 


and d/2 = {(|, |) : (u, u) S d}. Here, with a slight abuse of notation, we have 

Tn n (A\Ai/ 2 ) = {d : d S Tn and d C A\Ai/ 2 }- 

Note that all elements in T„ are shape regular, i.e. the radius of the largest 
inscribed circle of d is of the same order 0(k) as the diameter of d. Figure 
depicts the locally refined mesh for n = I, 2, 3. The meshes are defined in such 
a way, that the left half of the mesh T„ is a scaled version of the mesh Tn_i. 


This refinement scheme is presented in more detail in Takacs (2012). 





Figure 2: Meshes Ti, T 2 and T 3 


3.2 Hierarchical function space 

Now we can define a (hierarchical) mapped B-spline space over this mesh. The 
following definition is based on the construction presented in Takacs (2015). 
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Definition 2. The basis S„ over ]0, Ip is defined via 

§n = = b'^{s)Bj{t) : 1 <i <p+l and 1 < j < i + l| 

{Kj ^ b^{s)b]{t) ■. p+l + l<i<p + 2 and 1 < J < p + 2| 

“ b'^{s)b'j{t) :p + 2 + l<z<p + 4 and 1 < j < p + dj 

U Ihlj = b^{s)b'^{t) : p + 2"-i + 1 < i < p + 2" and 1 < j < p + 2”| , 

where Bj{t) is the j-th Bernstein polynomial of degree i. This defines a function 
space >V„ over A via 

Wn = span = S”^- o u“^ : G §„| . 

Lemma 1. The space Wn and the corresponding basis fulfill the following prop¬ 
erties. 

(A) The space Wn is a subspace ofWn 

(B) For all 6 G T„ the space fulfills C Wn p C Qp , where P^ is the space of 

bivariate polynomials of total degree < p and Qp is the space of bivariate 

polynomials of maximum degree < p. 

(C) The basis functions fulfill 

= fil-\2u,2v) 

for all 1 < i < 2"“^ and for all j. 

(D) The basis forms a partition of unity. 

We omit the simple proof of this lemma. 


3.3 Stable projection operator 


Let no = [log 2 ( 8 p)]. Then we have /iq = 1/2"“. In the following we consider 
n > uq. This is necessary for the proofs of Lemma and Theorem to be able 
to properly split the contributions from the singular point at the left and the 
regular part at the right of the parameter domain. We consider the standard 


dual basis for univariate B-splines, see e.g. Schumaker (2007); da Veiga et al. 
(2014) for more details. 


Definition 3. Let : L^(0,1) —>■ K be the dual basis for 6/ with A^(6”) = (5/, 
for z, fc G {1,..., 2" + p}. 
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Moreover, consider a dual basis for polynomials on Ah„. 
Definition 4. Let : L‘^{Ahg) —>■ M be dual functionals, with 


for fc = 1,... ,p + 1 and Z = 1,..., fc. 

We assume that all the functionals in Definitions [Hand [Hare bounded in L^. 
This is no restriction, see Schumaker ( |2007[ ). 

Let cl be the mapping from —>■ with (u,u) i—>■ {cu,cv). Then we can 

define a global dual basis for all n > no. 

Definition 5. Let n> no- Given tp G L^(A), we set 




for all fc = 1,... ,p + 1 and all Z = 1,..., fc. Moreover, we set 

G V o u) 

for all Zc = p + 2,..., 2” +p and all Z = 1,..., +p. Here m(k) = |’log 2 (Ze — 

p+ 1)1. 

It is easy to see, that ^ is in fact a dual basis for the basis on A. 

Definition 6. The dual functionals naturally define a projection operator from 
L^(A) onto >V„, with 

2"+p 

nvvjv^)= E E'^mW^ 

fc=i I 


where we appropriately sum over all Z. 

Before we prove our main Theorem [^ we recall two preliminary results. 

Lemma 2. Let n > ng. There exists a constant Cu > 0, depending only on p, 
such that for all (p G L'^(A\Afig) and for all S G T„ n (A\A 3 / 8 ) we have 




LHS) 


A: Cr ||y)|| ^2(5) 


( 1 ) 


where 8 is the support extension of 6 as presented in 


Bazilevs et al. 


(2006). 
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Proof. A proof of this lemma follows directly from the standard theory for 
regularly parameterized domains (see e.g. Bazilevs et al.|(2006|; da Veiga et al. 


(2014)), since A\Ahg is parameterized via Uq : ]ho,l[ x ]0,1[ — A\Aho and 
the support extension of any 6 G A\A 3 /g remains in A\A( 3 / 8 _p^). Since Uq 
is chosen large enough, we have ^ — ph > ho for all h < ho, and therefore 
S C A\Ahg. In that case, the constant Cr depends on the degree p and on the 
Jacobian determinant of the mapping Uq, which depends on no only. Since n-o 
is fully determined by p, the constant Cn depends on p only. □ 


On the refinement level uq it is obvious that the operator is bounded. 


Lemma 3. There exists a constant Co > 0, depending only on p, such that for 
all ip G L^(A) and for all S G we have 


n 




L2(5) 


^ C'o I|(/5||i2(,5) 


( 2 ) 


Proof. This lemma follows directly from the boundedness of the dual function¬ 
als. There is no condition on the scaling needed, since we consider the coarse 
level no only. The constant Oo depends on p and no only. Again, no is fully 
determined by p, which concludes the proof. □ 


Now we can prove the main result, the stability of the projection operator. 


Theorem 7. Let n > uq. There exists a constant Cs > 0, depending only on 
p, such that for all p G L^(A) and for all S G T„ we have 




(3) 


Proof. We split the proof into three parts. 


(A) Let S C A\A 3 / 8 . Then this statement follows directly from Lemma 

(B) Let S C A 3/8 and there exists an to € Z+, with m < n — no, such that 
2’"I((5) C A 3 / 4 \A 3 / 8 - Due to Lemma[^(C) we have 

bVniAa/s = bVn-l|A3/4 ° 21 


for all n > log 2 ( 8 p). Hence we conclude 
2 




L^S) 


= 2 “ 


n 


w„_. 


(<po(2-“I)) 


L2(2™I(5)) 


The desired result follows from Lemma |2] and because 






















(C) We have 




L^S) 


= 2’^°~ 


n 


w„ 


(^o(2"»-”I)) 


L2(2"-"oI(5)) 


Here, the bound is fulfilled because of Lemma 


This concludes the proof with Cs = max(Cij, Co). 


□ 


Having a stable projection operator, we can prove approximation error bounds 
on n for a special class of geometry mappings G. But before we present the 
proof, we discuss a possible generalization. 


3.4 Generalization to locally quasi-uniform knot vectors 

The results can be generalized to certain configurations of non-uniform refine¬ 
ments. Let T„ be a sequence of locally quasi-uniform meshes with Tq = {A}. 
In that case one can also show that there exists a bounded operator, project¬ 
ing onto the piecewise polynomial function space over the locally quasi-uniform 
mesh. The stability constant then depends on the degree p and on the quasi¬ 
uniformity constants, i.e. the maximum ratio between the diameter of an el¬ 
ement and the diameter of the inscribed circle, as well as the maximum ratio 
between the diameters of neighboring elements. Figure [^depicts a sequence of 
three locally quasi-uniform meshes. 




This generalization is based on the fact, that Lemma is actually valid for 
any mesh on A\A 3/8 (see da Veiga et al. (2014|). Moreover, Lemmacan be 
generalized to arbitrary initial meshes. Then the constant depends only on the 
degree and on the quasi-uniformity of the initial mesh. We do not go into the 
details of this generalization here, but continue with the proof of approximation 
properties in the uniform case. 
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4 Approximation error bounds 


We first prove approximation error bounds on the triangle A and then show 
bounds on ft, using the equivalence of norms presented in Assumption 


4.1 Approximation on A 

On A we have the following approximation error bound. 


Theorem 8. Let n > uq. Then there exists a constant C > 0, depending only 
on the degree p, such that for all (p G for all q < p + 1, and for all 

S G Tn we have 




Hi{S) 


< ChP-‘i+^\ip\ 


Hp+qi)- 


Proof. We show this result only for q = 0. The statement for general q follows 
immediately. According to Bramble and Hilbert (1970), there exists a constant 
Ca, depending only on p, and a polynomial gVp, such that 


The constant is uniformly bounded, due to the quasi-uniformity of the support 
extension S of the element <5. Due to Lemma[^(B) we have ip G Wn, hence 




LHS) 




2 

L^(S) 


From Theorem [3 it follows that 




Hence, we conclude 




LHS) 


< (1 + C|) \\ip — 1p\\ ^2(1) < (1 + C'g)C'Ahf^^‘^\ip\ 




which concludes the proof with C = Ca\/^ + Cg. 


□ 


We can now apply this theorem to prove bounds on the mapped physical 
domain H. 


4.2 Approximation on 12 

This bound can be extended directly to some geometry mapping G fulfilling 
additional regularity criteria. 
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Theorem 9. Let n > inax(no, n*), Zei G = F o u, where F = ;^(Fi, ^ 2 )^ with 
Fq,Fi,F 2 G W„» and det VF > c > 0, and let 


nv„ iv) = 


n^J(cpoF).Fo) 

Fn 


oF“ 


Then there exists a constant C > 0, depending only on the degree p and on the 
geometry parameterization G, such that for all (p G for all q < p + 1, 

and for all uj = F((5), with 6 G T^, we have 


W <ChF 

Proof. The proof of this theorem is a direct consequence of Theorem and the 
equivalence of norms as stated in Assumption The assumption is fulfilled 
because of Lemma (A) and Proposition □ 

Note that the condition on the geometry mapping G is a real restriction. A 
general NURBS mapping G cannot be represented as a composition of a singu¬ 
lar bilinear mapping u with a regular mapping F = -^{Fi, F 2 )'^ which fulfills 
Fo,Fi,F 2 G yVn- Hence, for general NURBS parameterizations, the equivalence 
of norms in Assumption is not fulfilled. However, numerical evidence shows, 
that in many cases the equivalence of norms is not necessary for optimal order 
of approximation. This phenomenon may be studied in more detail. 


5 Conclusion 

In this paper we could prove approximation error bounds for isogeometric func¬ 
tion spaces over singularly parameterized domains, thus extending known results 
in isogeometric analysis to singular patches. We show these bounds for a hierar¬ 
chically refined subspace of the full tensor product spline space. The refinement 
is considered to be uniform with h, which is not a necessary condition, as we 
pointed out in Section [3^ We however restricted our study to a certain type of 
singular patches that are derived from regularly mapped triangular domains. In 
the framework we presented, this is necessary in order to guarantee the equiv¬ 
alence of norms between the triangle A and the mapped domain U. It is not 
clear yet, how these results extend to more general conhgurations where this 
equivalence is not satisfied near the singularity. 
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